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Abstract. We present a study of 40 low-latitude unidentified 3EG gamma-ray 
sources which were found to be not positionally coincident with any known class of 
potential gamma-ray emitters in the Galaxy (Romero, Benaglia & Torres, 1999). 
We have performed a variability analysis which reveals that many of these 40 
sources are variable. These sources have, in addition, a steep mean value of the 
gamma-ray spectral index, < F 2.41 ± 0.2, which, combined with the high 
level of variability, seems to rule out a pulsar origin. The positional coincidences 
with uncatalogued candidates to supernova remnants were also studied. Only 7 
sources in the sample are spatially coincident with these candidates, a result that 
is shown to be consistent with the expected level of pure chance association. 
A complementary search for weak radio counterparts was also conducted and the 
results are presented as an extensive table containing all significant point-like radio 
sources within the 40 EGRET fields. We argue that in order to produce the high 
variability, steep gamma-ray spectra, and absence of strong radio counterparts 
observed in some of the gamma-ray sources of our sample, a new class of objects 
should be postulated, and we analyze a viable candidate. 
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1. Introduction 



The main results of the very successful Energetic Gamma Ray Telescope (EGRET) of 
the Compton Gamma Ray Observatory (CGRO) are contained in the Third EGRET 
(3EG) catalog (Hartman et al. 1999), which includes observations carried out between 
22 April 1991 and 3 October 1995, and lists 271 point sources. Of these sources, 170 
have no conclusive counterparts at lower frequencies. In the latitude range |&| < 10°, 
there are 81 unidentified sources, a number which doubles that reported in the previous 
2EG catalog and its supplement (Thompson et al. 1995, 1996), and clearly confirms the 
existence of a population of gamma-ray emitters in the Galaxy. 

The nature of this population of low-latitude gamma-ray sources is not clear. Several 
types of objects have been proposed as possible counterparts in the past: supernova 
remnants (e.g. Sturner & Dermer 1995, Esposito et al. 1996, Combi et al. 1998a), massive 
stars with strong winds (e.g. Montmerle 1979, Romero et al. 1999, Benaglia et al. 2000), 
young pulsars (e.g. Yadigaroglu & Romani 1997), and compact clouds in star- forming 
regions (e.g. Casse & Paul 1980). Gehrels et al. (2000) have recently shown that there 
is a population of sources associated with the Gould belt, as it was originally suggested 
by Grenier (1995). This suggests that young stellar objects play an important role in the 
generation of the observed emission. 

In a recent paper (Romero, Benaglia & Torres 1999, hereafter Paper I), some of us 
presented a study of the level of positional correlation between gamma-ray sources in the 
3EG catalog and a variety of galactic objects. A statistically very significant correlation 
was found in the case of SNRs and OB star associations. In that work it was also found 
that a group of 42 low- latitude 3EG sources do not present a likely counterpart. Two of 
these sources were finally discovered to be artifacts related to the high intensity of the 
Vela pulsar (Hartman 2000). Very recently, Zhang et al. (2000) proposed that most of 
the sources spatially coincident with SNRs and OB associations could be pulsars. They 
have presented a variability study of the 38 3EG sources that are coincident with objects 
of these classes, arguing that most of them are not variable. 

In this work we focus our attention on the group of 40 gamma-ray sources which 
are not coincident with any known galactic object thought to be capable of producing 
gamma-ray emission. Classical galactic counterparts are expected to be non-variable on 
timescales of a few years and, consequently, we shall discuss the time behavior of the 
sources in our sample. The structure of the paper is as follows. In Section 2 we define the 
sample. A variability analysis is presented in Section 3. Section 4 deals with the possible 
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positional correlation between sources in our sample and SNR candidates not included 
in Green's (1998) catalog. In Section 5 wc explore possible weak radio counterparts of 
the gamma-ray sources and present an extensive table with all significant point-like radio 
sources within the 40 EGRET positional error boxes that, we expect, will be very helpful 
for future research in this field. In Section 6 we discuss the case for isolated Kerr-Newman 
black holes as possible generators of some of the gamma-ray sources, and finally, in Section 
7, we draw our conclusions. Three appendices present detailed tables with data used in 
our research. 

2. Sample 

From the original set of 81 (75 if we exclude the six sources thought to be artifacts 
produced by the Vela pulsar) unidentified sources at |6| < 10°, an apparently large 
fraction of 50% (40 sources) remain without any known likely galactic counterpart. 
We list these sources in Table 1, where we give their 3EG name, their galactic coordinates, 
the error in their position -assumed to be represented by the 95% confidence contours 
given by the 3EG catalog-, and their gamma spectral index T (such that the photon 
distribution is given by N{E) oc E~^). We also present four further columns related to 
the variability analysis that will be explained below. A f-symbol on the source name 
indicates a source which was suggested as possible AGN by Hartman et al. (1999). It is 
worth noticing that the mean value of the spectral index is quite steep for these sources: 
2.41 ±0.2, steeper that the steepest pulsar spectrum known (Thompson et al. 1994). The 
distribution of the spectral index is shown in Figure 1. Only three sources out of 40 have 
r < 2. 

3. Gamma-ray variability 

We shall now assess the possible long term variability (on time-scale of years) of the 
sources in our sample. We will adopt the following criterion. We define a mean weighted 
value for the EGRET flux as: 



Here, A^vp is the number of viewing periods for each gamma-ray source. Wc take into 
account only single viewing periods (i.e. periods like 305+, virgo, or the combined P12, 
P123, and P1234 were not taken into account). F{i) is the observed flux in the i*''-pcriod, 
whereas e(i) is the corresponding error in the observed flux. We have taken these data 
directly from the 3EG catalog. For those observations in which the significance {\/TS in 
the EGRET catalog) is greater than 3(t, wc took the error as e{i) = F{i)/\/TS. However, 
many of the observations are in fact upper bounds on the fiux, with significance below 




n -1 
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3cr. For these ones, we assume both F{i) and e{i) as half the value of the upper bound. 
We then define the fluctuation index fj, as: 

/i = 100 X CTsd X {F)-^ . (2) 

In this expression, cTgd is the standard deviation of the flux measurements, taking into 
account the previous considerations. 

Results for the 40 sources are given in Table 1, where we show four columns corre- 
sponding to (F) , (Tsd, Nyp and /j,. 

In order to remove as far as possible spurious variability introduced by the observing 
system, we computed the fluctuation index /i for the confirmed gamma-ray pulsars in 
the 3EG catalog (see Cusumano et al. 2000 for the latest identification). We adopt the 
physical criterion that pulsars are -i.e. by definition- non-variable gamma-ray sources. 
Then, any non-null /U- value for pulsars is attributed to experimental uncertainty. We then 
define an averaged statistical index of variability, I, as 

J /^source Msource 

^ M ^pulsars 26.9 



Table 2 shows the variability results for the pulsars in the 3EG catalog, whereas the 
last column of Table 1 shows the results of the variability index / for our sample of 
unidentified 3EG sources. 

The adopted variability criterion is then to consider that any source for which the 
/i- value is less than the upper value of ji for pulsars (/imax = 40.5), is a non- variable 
source. At least, w(; can assure that for this kind of source, and with the data now at 
hand, it is not possible to discriminate any significant long-term variation in its gamma 
emission. We shall also consider that sources with /i-values between 40.5 and 65 are 
dubious cases, and that sources with /j, > 65 are variable. In terms of the averaged index 
/, this is equivalent to saying that variable sources will be those with / > 2.5, which is 
3cr away from the statistical variability of pulsars. The value of la is naively obtained 
from the standard deviation of the /-index for pulsars. Since this is a population of just 
a few members, this forces us to be careful in assigning to a source the "variable" status. 

We adopted this criterion, previously used in blazar variability analysis by some of 
us (Romero et al. 1994) and applied to 3EG sources by Zhang et al. (2000), rather than 
the similar one used by McLaughlin et al. (1996) and Wallace et al. (2000) in order to 
allow a direct comparison with the spurious statistical variability shown by pulsars. Since 
the /-index establishes how variable a source is with respect to the pulsar population, 
it can be considered as indicative of how reliable we can be about the possible physical 
variation in the gamma emission. The combination of the results of Zhang et al. (2000) 
with those presented in this paper yields the variability index for all unidentified low- 
latitude sources in the 3EG catalog, see Appendix C. Note, however, that Zhang et 
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al.'s study differs from tlie present one in that they considered PSR B1951-I-32 for the 
determination of < /x >puisars- This pulsar is not in the 3EG catalog. The fluxes used by 
Zhang et al. were taken from data reduced using a different technique. The inclusion of 
this pulsar makes < /i >puisars= 77.7 ± 50.0 instead of the much lower value we obtain. 
This high value of < /x > pulsars causes the /-index for AGNs to be compatible with a non- 
variable population, and it is also the reason why the majority of Zhang et al.'s sources 
are non-variable ones. In Appendix C, the values of I for all sources at low latitudes 
< 10") are computed under our normalization. 

A more comprehensive study on the long term variability (more than one month) 
of the gamma-ray sources was recently presented by Tompkins (1999). This work re- 
analyzed the EGRET data to take into account not only all sources included in the 3EG 
catalog, but also the 145 marginal sources that were detected but not included in the 
final official list. The maximum likelihood set of source fluxes was then computed. From 
those fluxes, a new statistics measuring the variability is defined as r = cr//x, where cr 
is the standard deviation of the fluxes and n their average value. The strength of this 
approach lies in that it takes into account possible fluctuations from the background and 
from neighboring sources, careful sensitivity corrections throughout EGRET lifetime, 
and others systematic errors. On the other hand, our method can be very useful as a first 
approach: when one takes a fairly safe assumption as the threshold for variability (such 
as the one introduced above, i.e. a source is classified as variable if / > 2.5) it yields 
results that are compatible with Tompkins'. For instance, if we consider the sources that 
Tompkins found as the most variable ones, we find that only two of them are included 
in our sample. These have / = 2.96 and I = 5.33, respectively. If we now consider the 
sources that Tompkins classified as less variable, we find that two of them are in our 
sample, with indices / = 1.00 and / = 1.86. The rest of the cases are compatible too, 
but loss strength should be put on the results, since most of them are dubious cases 
in both schemes. We then found that our method, avoiding most of the criticism of 
previous studies (Tompkins 1999) is well-fitted as a discriminator between variable and 
non-variable sources using the publicly available EGRET data. 

With our adopted criterion, we find that only a low 25.0% of the sources in our sample 
are non-variable, 42.5% are dubious cases, and 32.5% are variable. In this latter group, 
three sources have a variability index 5(t away from the statistical variations of pulsars. 
In Figure 2 we present the histogram for the variability index / in our sample, along 
with a Gaussian fit which presents its peak at / =2.0 and a standard deviation equal to 
(7/ =0.7. In order to better evaluate these results, we have also computed the variability 
index for the 66 confirmed AGNs in the 3EG catalog. A table with these results (with 
the same columns as in Table 1) is given in Appendix B. Fifty-eight AGNs appear to be 
variable sources. We show the distribution of the variability index for AGNs in Figure 3. 
The Gaussian fit in this case peaks at / =2.3 and the standard deviation is cr/ =0.8. 
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Since the distribution of the unidentified sources under study is similar to that obtained 
for AGNs, which is a well-known variable population, we can conclude that many of our 
EGRET sources may actually vary on time-scales of years. 

In Fig. 4 we plot the computed value of / versus the observed value of the spectral 
index F, for the sources in our sample. We also show in the same figure the / - F 
relation for AGNs. Two solid horizontal lines represent the lower and upper bound for 
the statistical variability of the pulsar population. An horizontal dashed line stands for 
the value of the variability indices above which the variability of the sources can be 
established (i.e. above Scr level). From these plots, it appears that more than a single 
population is present. Pulsars should be located at the lower left corner of the figure: 
they have a hard spectral index and are non-variable. A different group of sources, with 
spectral index in the range ^ 2.1— 2.6 and being marginally variable or variable, 
appears towards the center of the figure. Most of the AGNs are located in this region of 
the plot. Finally, some sources seem to show a steepening in their spectral index with 
increasing values of variability. A Spearman correlation test for the variable sources yields 
a probability of about 8% for this being a pure chance effect. These sources could then 
be members of a new population of galactic objects or, perhaps, they might be AGNs 
seen through the disk of the Galaxy. 

Since AGNs are isotropically distributed, we can determine a mean value for the num- 
ber of AGNs detected by EGRET per square degree: TVagn = 8.9 x 10"'' AGN/degree^. 
At latitudes in the range —10° < b < 10°, the number of unidentified sources 
that are marginally variable or variable yields a number density Nunid = 4.3 x 10^'^ 
objects/degree^. This shows that there are five times more variable objects within the 
galactic plane than away from it. Then, we have to conclude that if the distribution of 
AGNs is isotropic, we have an excess of variable galactic sources representing a new class 
of population, as it was already noticed in previous studies by McLaughlin et al. (1996), 
by Mukherjee et al. (1997), by Tavani et al. (1997), and very recently by Wallace et al. 
(2000). Q 

Considering the sample of the 2EG catalog, Ozel and Thompson (1996) argued that 
one half or more of the unidentified EGRET sources lying at high latitudes are of galactic 
origin. Nice and Sawyer (1997) found no evidence that pulsars were the constituents of 
this putative galactic population. This agrees with our results, since we show here that 
this new class of galactic objects should resemble, from the point of view of gamma-ray 
observations, the AGN population. 

It is perhaps appropriate to come back to the discussion of what is a real source 
(as opposed to possible noise variations) in this study. We are considering that the 3EG 
catalog is composed of real sources, disregarding the level of flux detected in each case. 

^ The latter work focuses only on very short timescales. 
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The exception to this rule is the six artifacts related to the Vela pulsar, which disappear 
in a map where the Vela emission is suppressed (Hartman et al. 1999). However, since 
many detections in the EGRET catalog are below the 3a significance, and the fluxes 
are compatible with zero, a more careful treatment is in order. We can consider only 
those sources with ^TS > ber, which implies a flux above 30 xlO~* ph cm~^ s~^. In 
that situation we are left only with 16 sources, four of them having / > 2.5. One of 
these sources is 3EG 1828+0142, with / = 5.33, to which wc devoted a separate work 
elsewhere (Punsly et al. 2000). If the flux threshold is lowered to 18 x 10~^ ph cm~^ s~^, 
we are left with 35 sources, and now 11 of them have / > 2.5. Within this group is 3EG 
1735-1500, which has / = 8.86. It is then apparent that variable sources exist even when 
the flux threshold is raised. 

The relevance of the previous point lies mainly in the galactic longitude distribution 
of the 40 3EG sources of Table 1. In the left panel of Fig. 5 we show the histogram 
distribution with galactic longitude of all sources in our sample; it is clear that it is 
peaked towards the galactic center and the galactic anti-center. Although it is true that 
the EGRET exposure was greater in these regions, it appears that the number of sources 
of this sample that were detected at these positions is relatively larger than what results 
when the whole 3EG catalog is considered. If we take a more strict limit on the flux, 
such as given by taking sources only with ^/TS > 5a, we obtain the distribution shown 
on the right panel. That is, a more restrictive limit on the adopted minimum flux will 
preferentially delete the anti-center sources, leaving strong detections towards the galactic 
center, which is certainly more compatible with a real galactic population. 

It seems to be possible that at least a part of this excess in the number of sources 
towards the galactic center could be the low-latitude tail of the new population found by 
Gehrels et al. (2000) at mid-latitudes. This speculation is supported by the fact that the 
average spectral index of the sources associated with the Gould belt is 2.49 ± 0.04, quite 
compatible with the indices of the sources in the sample here considered. However, the 
sources in the mid-latitude population are weaker than many of the variable and strong 
ones of our group. If stars are responsible for the weaker sources, as suggested by Gehrels 
et al. (2000) and by Benaglia et al. (2000), then a different type of object, intrinsically 
much more luminous at gamma-rays, should be behind the stronger low-latitudes EGRET 
detections. 

A preliminary conclusion of this section could then be stated as follows: A new type 
of galactic objects may be needed to explain the behaviour of some unidentified 3EG 
sources at low latitudes. This population should be able to display high luminosities, 
steep gamma-ray spectral indices and significant gamina-ray variability. Fig. 4 suggests 
that our sample, although mostly of galactic origin, mimics some observational properties 
of AGNs. 
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4. Uncatalogued supernova remnants 

We consider in this section whether some of the sources in our sample may be associ- 
ated with recently proposed candidates to supernova remnants, not catalogued by Green 
(1998). Our interest in this search resides in the fact that young stellar objects like re- 
cently formed black holes and pulsars can still be associated with the gaseous remnant 
of the original supernova that created them. 

The diffuse non-thermal emission of the galactic disk, formed by the interaction of the 
leptonic component of the cosmic rays with the galactic magnetic field, is surely veiling 
many remnants of low surface brightness. Recent observational studies using filtering 
techniques in the analysis of radio data have revealed several new SNR candidates that 
are not yet included in the latest issue of Green's catalog, which was used in the study 
presented in Paper I (e.g. Duncan et al. 1995, Combi & Romero 1998, Combi et al. 1998b, 
1999a, Jonas 1999). In general, these new candidates are much more extended than those 
previously known. There are 101 of these weak non- thermal structures detected so far 
in the Galaxy. This number significantly extends Green's (1998) catalogue. The list of 
these new candidates, compiled from the papers and studies above mentioned is being 
published electronically as an Appendix. 

With the aim of finding the positional coincidences between our sample of 40 3EG 
unidentified sources and the uncatalogued candidates to SNRs, we used the code de- 
veloped in Paper I. We have found that only 7 gamma-ray sources in our sample are 
positionally coincident with non-thermal radio structures. The positional coincidences 
thus obtained are shown in Table 3, where we give the name of the gamma-ray sources, 
the central position of the SNR candidates, the distance from the EGRET source to the 
centre of the radio structure {6), and the addition of the radius of the extended sources 
and the uncertainty in the position of the EGRET detections (A). 

To estimate the statistical significance of these coincidences, we have numerically sim- 
ulated a large number of synthetic sets of EGRET sources using the code described in 
Paper I. As in that paper, the simulations were constrained to preserve the original gra- 
dient in the number distribution of unidentified gamma-ray sources towards the galactic 
plane. The results of this study show that the expected number of chance associations is 
10.4 ± 2.7, quite compatible with, and even bigger than, the actual result. We conclude, 
then, that there is no statistical evidence suggesting that the BEG sources analyzed in 
our sample are associated with known or potential SNRs. 

5. Search for extended radio counterparts in selected sources 

We have especially examined the radio sky towards the best estimate positions of the 
unidentified 3EG sources J0241+6103, J0435+6137, J0628+1847, J1014-5705, J1631- 
4033, J1735-1500, J1828+0142, J1928+1733, J2035+4441 and J2206+66. From Table 1, 
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it can be seen that these sources can be grouped into two sets, at the extremes of the 
variabihty and spectral index ranges. On one hand, we have six sources with low variabil- 
ity and relatively hard spectral indices. These sources could, in principle, be associated 
with pulsars or with unknown new SNR candidates. Only one of them, 3EG J1631-4033, 
presented positional superposition with a known candidate to SNR (see Table 3). There 
are, on the other hand, four sources with high variability (J-values ranging from 3.3 to 
8.8) and steep spectral index (F ranging from 2.08 to 3.24). 

We used continuum radio data at 1420 and 2400 MHz, from the surveys by Reich 
& Reich (1986), and Duncan et al. (1995), in order to study the possible presence of 
weak and extended radio sources within the EGRET positional error boxes of this set. 
Specifically, we searched for yet undetected low-brightness SNRs that could be associated 
with the pulsars or the black holes created in the original supernova explosion. We remove 
the background diffuse radiation using the Gaussian filtering technique described in a 
recent series of papers by Combi and coworkers (e.g. Combi & Romero 1995, 1998; 
Combi et al. 1998a, b, 1999a,b). The reader is referred to these papers for details. Three 
of the fields that were studied presented features interesting enough to deserve individual 
description. Below we briefly comment on two of them: 3EG J0435-I-6137 and 3EG J1631- 
40. The source 3EG J1828-I-01 is discussed separately in the paper by Punsly et al. (2000). 
The remaining fields were empty of extended radio structures or too confused even after 
background subtraction as to provide reliable information. 

5.1. 3EG J0435+6137 

In Fig. 5, upper panel, we show the filtered radio image of the region around 3EG 
J0435-I-6137 at 1.42 GHz (with an rms noise of 20 mK and angular resolution of 
34') super-imposed to the EGRET probability contours. Two apparently extended ra- 
dio sources are visible in the field. The strongest one, 87GB J0443-I-6118, is par- 
tially outside the 95% confidence contour of the gamma-ray detection. Another source, 
named 87GB J0441-I-6145, is located 20' away from the best estimated position of 
3EG J0435-I-6137, well within the 95% confidence contour. 

In Fig. 5, lower panel, wc show the central region of the radio field at 5 GHz and a 
resolution of 3.5' (Condon et al. 1989). Several point-sources arc visible in the frame. The 
identification of the main ones is also indicated in the figure. None of these sources is 
known to be cxtragalactic. One of them, 87GB J0442+614G, is coincident with the X-ray 
source IRXS J044239. 3+61404. Another source, 87GB J0435+6137 (also known as TXS 
0431+615), is located exactly at the best estimated position of 3EG J0435+6137. This 
source has a steep spectral index of ^ —1.2. 

Wc have made Gaussian convolutions of those groups of point sources that are nearly 
coincident with the apparent extended sources seen in the lower resolution map of the 
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upper panel, using a beamwidth of 34'. We have found that both apparent structures 
are likely artifacts resulting from the merging of the several point-like sources in the 
region and that, consequently, no extended non-thermal source can be associated with 
3EG J0435-I-6137. The counterpart, if it exists at all, should be one of the point sources. 
The most promising candidates seems to be 87GB J0442-I-6140, which is also detected 
at X-rays, and 87GB J0435-I-6137, which is the strongest one at 5 GHz. 

5.2. 3EG J1631-4033 

Figure 6 shows the filtered radio continuum map of the field around 3EG J1631-4033 
at 2.4 GHz (with an rms noise of 12 mJy bearn"^ and an angular resolution of 10.4'). 
The confidence contours of the likelihood test statistics of the EGRET detection are also 
shown as a superposed gray-scale. Four point-like radio sources can be identified with 
PMN sources (Griffith & Wright 1993), namely (listed by increasing Galactic longitude), 
PMN J1627-3952 (noted as SI in Figure 6), PMN J1631-4015 (S2), PMN J1636-4101 (S3), 
and PMN J1631-3956 (S4) . The angular separation between the location of the highest 
likelihood test statistic for the gamma-ray source and the position of PMN J1631-4015 
is just 22 arcminutes. The strongest source, located at {l,b) ~ (341.9°, -1-4.15°), can 
be identified with PMN J1636-4101. We summarize the radio information on the point 
sources in Table 4. Any of these sources can be considered as a potential counterpart. 
Additional and weaker sources in this region are listed in Table 5 (see below). 

Another interesting feature revealed by the radio map is the existence of a minimum 
in the continuum emission towards the centre of the field at {I, b) ^ (341.2°, -|-4.6°). This 
might be the relic of an old explosive event in the ISM in this region. 

6. Point-like counterparts: systematic search 

We have searched CATS (Vcrkhodanov ct al. 1997), SIMBAD and NED databases for any 
possible radio counterparts for the gamma-ray sources in our sample. We have prepared 
a table containing the results of this search. This table (Table 5, which is published only 
in electronic form due to its length) lists all significant point-like radio sources, with 
flux densities above 30 mJy at least at one wavelength, found within the 95 % EGRET 
confidence contours of the 40 gamma-ray sources of Table 1. 

In Table 5 there arc entries for each gamma-ray source listing, from left to right, 
radio survey used, source name, right ascension and declination with the corresponding 
uncertainties, observing frequency, flux density and its error, and comments. The different 
flelds are separated, and within each field (listed in increasing observing frequency), 
a different paragraph corresponds to each radio source. All listed sources (about 600, 
ranging from just a few to almost 50 per 3EG field) are located within the EGRET 95% 
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confidence contours. The material contained in this table could be very useful for future 
research in this field. 

7. Discussion: The black hole hypothesis 

The existence of a number of variable gamma-ray sources with steep spectra and without 
strong radio or X-ray counterparts can be explained postulating a galactic population 
of magnetized black holes (Punsly 1998a, b, 1999). In these objects, strong magnetic 
fields (~ 10^^ G) are supported by an orbiting charged ring or disk. The black hole itself 
is also charged and rotating (Kerr-Newman black hole), but the entire configuration 
has zero net charge, in such a way that it is not quickly neutralized by the accretion 
of difi^use interstellar matter. The resulting object has a magnetosphere similar to that 
presented by pulsars. However, because of the absence of a solid surface, thermal X- 
rays are not produced, as in a neutron star. In addition, there is no accretion disk 
capable of generating strong X-ray emission as is expected in microquasars (Mirabel 
and Rodriguez 1998). The radiation of the system, which is non-pulsating due to the 
alignment of the rotation and magnetic axes, is produced by synchroton and inverse 
Compton mechanism in two relativistic electron-positron jets which propagate along 
the rotation axis in opposite directions, as it occurs in AGNs. Contrary to the case of 
most AGNs, no redshifted emission lines should be present and, in contrast to BL Lacs, 
no weak nebulosity is expected in the optical band. Pmisly et al. (2000) have recently 
shown that, when pair annihilation effects are taken into account in the calculation of 
the spectral energy distribution presented by magnetized black holes, a steep gamma-ray 
spectrum is produced. Variability naturally results in the model as the effect of shocks 
and firehose instabilities which induce large changes in the Doppler enhancement factor 
of the relativistic jets (Punsly et al. 2000). 

At high energies, self-Compton losses provide most of the gamma-ray luminosity, 
reaching values in the range lO''^ — 10'^^ erg s^^ for black holes of a few solar masses 
and polar magnetic fields of ^ 10^^ G (see Punsly et al. 2000 for details). If such an 
object is relatively close (e.g. ^ 1 kpc), it could appear as a typical imidentified EGRET 
source with F ^ 2.5 or steeper. Although the spectral energy distribution is broad, only 
weak emission at longer wavelengths is expected. Radio counterparts, for instance, would 
present fiux densities on the order of a few to a few tens of mjy at 5 GHz. Interestingly, 
one would also expect a mild correlation in variable gamma-ray sources between the 
variability and spectral indices, as seemingly observed in Figure 4. See the work by 
Punsly et al. (2000) for details on the model. 

Apart from 3EG J1828+0142, the source 3EG J0241+6103 (2EG .10241+6119) was 
already suggested as a possible identification of a magnetized black hole by one of us 
(Punsly 1999a, b). The source 2EG J0241-I-6119 is known to be variable since the work 
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by Tavani et al. (1998), who discarded a pulsar origin. This source appeared in a search for 
gamma-ray emitters, with strong emission also in X-rays, but weak output in radio. The 
binary system LSI -|-61°303, which is within the 95% probability contour of the gamma- 
ray detection, contains an unseen compact object suggested as the possible generator 
of the observed gamma rays (Punsly 1999a, b). However, a positive identification with 
the 3EG detection is not conclusive because of the fact that there is no signature in 
the gamma emission that can be correlated with the binary dynamics. Note that the 
magnetized black hole model in Punsly (1999b) does not produce any modulation of the 
gamma-ray signal on orbital time scales to a first order. 

8. Conclusions 

In this paper we have shown that there exists a group of unidentified, low-latitude, 
gamma-ray sources with the following characteristics: 

1. They do not positionally coincide with any known galactic object such as Wolf-Rayet 
or Of massive stars, catalogued SNRs, and OB associations. The probability that a set 
of 40 sources, with the same positional errors as the EGRET ones, lacks any positional 
correlation with these galactic objects totally by chance is less than 3 x 10~^. 

2. Even considering positional coincidences with 101 identified candidates to SNRs, in 
addition to Green's (1998) catalog, only 7 out of 40 sources present correlation, which 
is consistent with being the mere effect of chance. 

3. The sample presents interesting variability features: 32 % are variable and 45 % of 
them are dubious cases. 

4. The sample presents a steep average spectral index < F >= 2.41 ± 0.2; only 3 out of 
40 sources have F < 2. 

5. Due to isotropy constraints on the average density of active galactic nuclei, not all 
of the 40 sources considered here can be AGNs seen through the disk of galaxy. The 
surface density of variable unidentified sources at low galactic latitudes is five times 
higher than that of AGNs out of the plane. This, and the fact that there are 5 AGNs 
already detected at |6| < 10°, makes that if the distribution of AGNs is isotropic, just a 
few out of 40 gamma-ray sources could be AGNs. Indeed, if we extrapolate the AGN 
number density out of the plane to latitudes within (—10°, 10°), we would expect 
6.5 AGNs to be detected by EGRET, leaving room for just a couple of remaining 
unidentified objects to be AGNs, one, perhaps, being B2013-I-370 (Mukherjee et al. 
2000). 

Altogether, these features arc not in accordance with usual models of gamma-ray 
emitters: they cannot be explained by pulsars, SNRs in interaction with nearby clouds, 
massive stars, or AGNs. In addition, wc have found an apparent trend for sources with 
a higher degree of variability to present steeper photon spectral indices. This latter 
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result, however, needs confirmation with the more complete sample to be obtained by the 
next generation of gamma-ray satellites (Gehrels et al. 1999). With the currently large 
positional error boxes, we have no hope to completely identify the source for all these 
detections, but if, after the launch of future gamma-ray satellites the previously quoted 
characteristics are confirmed, a new population of objects probably will be needed. 

To our knowledge, the only model capable of reproducing all the peculiar features 
presented by some of the sources in our sample (steep spectral gamma-ray index, high 
variability, absence of clear lower-frequency counterparts, etc) is that of isolated Kerr- 
Newman black holes (Punsly et al. 2000). This does not exclude the possibility that 
alternative models based on, for instance, stellar high-energy variability (Benaglia et al. 
2000), or pulsar abnormal activity (e.g. Tavani and Arons 1997), could play an important 
role explaining some detections. 

The hypothesis of a population of magnetized black holes in the Galaxy can be tested 
with the forthcoming GLAST and INTEGRAL satellites, which will allow investigators 
to increase the sample and measure the spectral energy distribution for particular sources 
with high precision. In the case of INTEGRAL, the spectrometer SPI and the imager 
IBIS could detect the electron-positron annihilation features of the black hole jets, as 
indicated by Punsly et al. (2000). 
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Appendix A: SNR candidates 

In this appendix we present the complete list of 101 candidates to SNRs that were used 
in our study of positional coincidences. This list is presented electronically in Table 6 and 
it consists of three columns: galactic coordinates (/, b) and the diameter of the extended 
source (9) in arcminutes. The positions of these sources were reported by Duncan et 
al. (1995), Combi & Romero (1998), Combi et al. (1998b, 1999a), and Jonas (1999). In 
general, these new candidates are much more extended and weaker than those previously 
known. 



Appendix B: Variability of AGNs 

In this appendix we present data for the variability indices of the 66 AGNs already 
detected in the 3EG catalog. These data were used both in Figure 4 and for making a 
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comparison with our sample in the main text. We provide electronicaUy (Table 7) the 
same columns as those given in Table 1 for the unidentified sources. 



Appendix C: Variability analysis of all low latitude gamma-ray sources 

If we combine the present work with that of Zhang et al. (2000), (taking into account 
the same normalization for < /i >puisars, as explained above) who studied the variability 
of the sources which were positionally coincident with OB associations and SNRs as 
discovered by Romero et al. (1999), we can now make an assessment of the variability 
of the whole sample of unidentified sources at low galactic latitudes. There are 81 such 
sources, six of them (3EG J0824-4610, 0827-4247, 0828-4954, 0841-4356, 0859-4257 
and 0848—4429) were recently reported to be artifacts produced by the high intensity 
emission of the nearby Vela pulsar. Another one, 3EG J0747— 3412, was found to be 
coincident only with a WR star by Romero et al. (1999), and its variability index is here 
reported for the first time to be / = 2.81. For future reference we present electronically in 
Table 8 the 3EG source name and the variability index / for all low-latitude unidentified 
sources, normalized to the gamma-ray pulsars detected in the 3EG catalog. 
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Fig. 1. Distribution of the spectral indices of the 40 unidentified sources. 
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Table 1. Unidentified 3EG sources without known possible galactic counterparts. 
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This source also displays short-term variability (Wallace et al. 2000). 
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Table 2. Statistical variability of pulsars. The last one was identified by Cusumano et 
al. (2000). 
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Table 3. Positional coincidences with candidates to SNRs. 
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Table 4. Characteristics of the point-like radio sources inside the 7-ray contour of 3EG 
J1631-4033 
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Fig. 2. Variability index distribution. The spline curve is the best Gaussian fit, it peaks 
at I =2.0 and presents a standard deviation of ai =0.7. These values suggest that many 
of the unidentified sources are variable. 
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Fig. 3. Variability index distribution for the 66 AGNs in the 3EG catalog. The spline 
curve is the best Gaussian fit, it peaks at / = 2.3 and has a standard deviation of crj =0.8. 
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Fig. 4. Variability versus photon spectral index plot. We also present the values for 
pulsars and 66 AGNs in the 3EG catalog. Horizontal solid lines represent the lowest and 
maximum values for pulsars. The dotted line corresponds to the limit above which all 
sources are variable {3a). 
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Fig. 5. Galactic longitude distribution of the 3EG sources considered in our sample, for 
different flux lower limits. See text for explanation. 
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Fig. 6. Upper panel: Map of 

the radio field around 3EG 
J0435+6f37 at f.42 GHz (an- 
gular resolution of 34') with 
the diffuse emission filtered us- 
ing a beamwidth of 120' x 
120'. Contours are labeled in 
steps of 0.02 K from 0.05 to 
0.22 K in brightness temper- 
ature. The superimposed grey- 
scaled levels represent the 99%, 
95%, 68%, and 50% statistical 
probability that the gamma-ray 
source lies within each contour. 
Lower panel: Image at 5 GHz 
of the central field obtained with 
the former NRAO 91-m telescope 
at Green Bank (Condon et al. 
1994). Contours are shown start- 
ing at 10 mjy beam^^, in steps 
of 10 mJy beam^^. 
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Fig. 7. Radio image of the field 
of 3EG J1631-40 at 2.4 GHz 
(angular resolution of 10.4') af- 
ter the subtraction of the diffuse 
emission. Superimposed on the 
radio image we show the confi- 
dence contours of the likelihood 
test statistics of the EGRET 
detection. The best result for 
revealing a counterpart of the 
gamma-ray source is obtained for 
a filtering beamwidth of 90' x 90' 
and 6 iterations. Radio contours 
are shown at steps of 0.05 from 
0.05 to 6 Jy beam^^. 



